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Abstract 

We consider statistics of the passive scalar on distances much larger than 
the pumping scale. Such statistics is determined by statistics of Lagrangian 
contraction that is by probabilities of initially distant fluid particles to come 
close. At the Batchelor limit of spatially smooth velocity, the breakdown of 
scale invariance is established for scalar statistics. 
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INTRODUCTION 



If an external pumping supplies the fluctuations of the scalar at some scale 
L, then the advection by a spatially inhomogeneous velocity field produces 
scalar fluctuations at all scales, both larger and smaller than L. In an incom- 
pressible velocity field, the flux of the scalar variance flows downscales, this 
direct cascade is quite well understood by now0^. From a general physical 
viewpoint, it is also of interest to understand the properties of turbulence at 
scales larger than the pumping scale. If only direct cascade exists, one may 
expect equilibrium equipartition at large scales with the effective tempera- 
ture determined by small-scale turbulence^. The peculiarity of our problem 
is that we consider scalar fluctuations at the scales that are larger than the 
scale of excitation yet smaller than the correlation scale of the velocity field, 
which provides for mixing of the scalar. Although we find simultaneous cor- 
relation functions of different orders, it is yet unclear if such a statistics can 
be described by any thermodynamics-like variational principle.. 

Since we are interested in the behavior of the passive scalar on large scales, 
the diffusivity can be neglected, so that the properties of the scalar statistics 
are solely due to Lagrangian dynamics. In a turbulent flow, the distances 
between fluid particles generally grow with time. The law of such growth 
determines the correlation functions of the scalar at the distances smaller 
than L. For example, the pair correlation function {6{r)9{0)) is proportional 
to the average time two fluid particles spend within the pumping correla- 
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tion scale. For r < that is the time when separation grows from r to 
L. On the contrary, the scalar statistics at scales larger than L is related 
to the probabilities of initially distant particles to come close. Study of the 
large-scale statistics thus reveals new information on the properties of La- 
grangian dynamics in a random flow. We shall show below that the statistics 
of Lagrangian contraction critically depends on the spatial smoothness of the 
velocity field. We shall argue that non-smooth velocity provides for a scale- 
invariant statistics of a scalar which is even getting Gaussian at the limiting 
case of extremely irregular velocity. Contrary, the statistics is rather peculiar 
at spatially smooth random flow (so-called Batchelor limit): it demonstrates 
strong intermittency and non-Gaussianity at large scales. Another unex- 
pected feature of the scalar statistics in this limit is a total breakdown of 
scale invariance: not only the scaling exponents are anomalous (i.e. don't 
grow linearly with the order of correlation function) yet even any given corre- 
lation function is not generally scale invariant (that is the scaling exponents 
depend on the angles between the vectors connecting the points) . 

The paper is organized as follows. We introduce the problem and discuss 
the results that could be understood qualitatively in Sect. |. These results are 
supported by straightforward calculations within the framework of Kraichnan 
modelH, presented in Sects. y-fV|. We briefly describe the case of non- 
smooth velocity in Sect. [V]. We consider arbitrary space dimensionahty d. 
Two-dimensional case deserves separate consideration due to an additional 
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degeneracy. 



I. QUALITATIVE DESCRIPTION 

The evolution of the passive scalar 9{r,t) under the action of velocity 
v{r^t) and pumping (J){r^t) is described by the equation 

dte + vVe = (j). (1.1) 

Let us introduce Lagrangian trajectories ^(r, t) determined by the equation 
dtQ = v{t, q) and by the initial condition Q{0,r) = r. Next, introducing 
9{t, r) = 9{t, q) we rewrite ([TTT]) as dt9 = 0, which gives the formal solution 

e{0,r)= f dt(t){t,Q). (1.2) 

Here we have taken into account that at i = the functions 9 and 9 coincide. 

Both V and (f) are assumed to be random functions of time and 
space. We will examine n-point correlation functions of the passive scalar 
F„(ri, . . . ,r„) = {9{ri) . . .9{rn)), averaged over both the statistics of the 
advecting velocity v and of the pumping 0. Since our main interest here is 
to study the scalar statistics on large distances and time scales, then without 
lost of generality we may consider pumping statistics to be close to white 
Gaussian 

(0(ti, n)0(t2, r^)) = 5{h - h)x (ki -r2\). (1.3) 

Here x is assumed to decay on a scale L. One can treat a deviation from 
Gaussianity by introducing the three-point pumping correlation function 
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{(J){ti,r,)(j){t2,r2)^{h,r^)) (1.4) 
= 5{ti - t2)5{ti - U)xi (|ri - r2|, \ri - rs], |r2 - rsl) , 

where X3 is supposed to have the same characteristic length L as x- Note, 
that even when xa introduces a small correction to the Gaussian statistics of 
the source, it produces a new effect, making the odd correlation functions of 
the scalar non-zero. The correlation functions can be represented as 



F2n= J dti... J dtn ixlRuiti)] . . . x[R2n-l,2n{tn)]) + • • • ' (l-^) 

— oo —oo 


F2n+1= J dh... J dtn{x3[Ruiti),Ri3{ti),R23{tl)] 

—oo ~oo 

XX[R45{t2)] . . . x[R2n,2n+l{tn)]) + ■■■ ^ (1-6) 

where angular brackets mean averaging over the statistics of the velocity 
and one should perform summation over all sets of the pairs of the points 
Vi. Using we have written the correlation functions in terms of the 

Lagrangian separations 

R,j{t) = \Q{t,r,)- Q{t,r,)\. (1.7) 

Most of this paper is concerned with the case where the velocity field can be 
considered spatially smooth, which means we can write 

Va{t,ri) - Va{t,r2) = (Jap{t)[rif3 - r2/?] . (1.8) 

Here cTq,^ is the random strain matrix depending only on time. At such a 
velocity field, the distances Rij{t) grow exponentially, the stretching rate 
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X{t) = ln[R{t) / R{0)] has Gaussian statistics with nonzero mean A and with 
the dispersion decreasing as ^^^/^ at time intervals far exceeding the correla- 
tion time of (7. 

Let us briefly remind the properties of the small-scale scalar statistics as 
they follow from (|i.b| - |i.8|) . When r <C the pair correlation function is pro- 
portional to the mean time when R{t) < L so that F2{r) = x(0)A^Mn(L/r) 
with logarithmic accuracyS'Si. With the same accuracy, the moments with 
n <C ln(L/r) are Gaussian at small scales^J^. 

The situation is drastically different at r > L. Now, non-zero correlation 
at two distant points appears only when two fluid particles manage to come 
there that were in the past within the pumping correlation length. We thus 
have to estimate the probability for the vector R{t) that was once within 
the pumping correlation length L to come exactly to the prescribed point 
r which is far away. Since the volume is conserved, then all the particles 
from the pumping volume will evolve in such a way to be stretched in a 
narrow strip with the length r. Assuming ergodicity [which requires that the 
stretching time A^^ln(r/L) is much larger than the strain correlation time] 
we thus come to the conclusion that the probability that two points separated 
by r belong to a "piece" of scalar originated from within L is given by the 
volume fraction {L/rY. That gives the law of the decrease of the two-point 
scalar correlation: F2 oc r~'^. 

The peculiarity of the smooth velocity field (|TT^) is that it preserves 
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straight lines under advection. That makes it easy to determine r-dependence 
of the correlation function of arbitrary order if all the points lie on a line. In 
this case, the history of stretching is the same for all the distances. Look- 
ing backward in time we may say that when the largest distance between 
points was within L then all other distances were as well. Therefore, the n- 
point correlation function for collinear geometry is determined by the largest 
distance: F„ oc r~'^. This is true also when different pairs of points lie on 
parallel lines. Note that the exponent is n-independent which corresponds to 
a strong intermittency and an extreme anomalous scaling. The fact that for 
collinear geometry ^ is due to strong correlation of the points along 
the line. 

When we consider an arbitrary geometry, the opposite takes place, namely 
the stretching of different non-collinear vectors is generally anti-correlated 
because of incompressibility and volume conservation. Indeed, for a smooth 
velocity field there exists a number of invariants, preserved by the flow. A d- 
volume eaia2...adPT ■ ■ ■ pT is conserved for any d Lagrangian trajectories Pi{t). 
In particular, for d = 2 there are area conservation laws eaf}PiP2 for any two 
vectors relating three points. Let us now consider a two-dimensional flow 
where the anti-correlation due to area conservation can be easily understood 
and the scaling for non-collinear geometry can be readily appreciated. Since 
the area of any triangle is conserved then three points that form a triangle 
with the area s much larger than LP' will never come within the pumping 
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correlation length. Therefore, the triple correlation function 



i^3(ri2,ri3,r23)= / dt (xsiRuih) , Ruih) , Risih)]) (1.9) 

— oo 

is determined by the asymptotic behavior of xs at ^ L that is very small. 

For example, if xs decays exponentially then F3 oc exp(— s/L^). On the other 

hand, for collinear geometry F3 oc r~^. We thus see that F3 as a function of 

the angle 1} between the vectors ru and ri3 has a sharp maximum at zero 

angle and decreases within the interval 'd ~ L^/r^ <C 1. 

Similar considerations apply for the fourth-order correlation function 


F,= J dh J dt2(x[i?12(tl)]x[i?34(t2)]) + ... , (1.10) 

—00 —00 

where dots stay for all possible permutations of points. Let us consider 
the contribution from the first term. Again, since the area \Ri2 x 7^34! is 
conserved, the answer is crucially dependent on the relation between \rs4 x 
rul and L^. When |r34 xru] <^ L'^ we have a collinear answer F4 oc r~^. Let 
us now consider the case of non-collinear geometry and find the probability 
of an event that during evolution R12 became of the order L, and then, at 
some other moment of time, reached L (only such events will contribute 
into F4). Note that, unlike the case of the third-order function, now there is 
a reducible part in pumping, which makes F4 nonzero (decaying as power of 
Vij) even when |r34 x ru] ^ L"^. The probability that Ru came to L is L'^ /r\2- 
Due to area conservation, there is an anti-correlation between R12 and i?34: 
if R12 ~ I/, than i?34 ~ ri2r^4/L. So probability for R^4 to come back to L is 
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L? / (ri2r34/L)^ = L^/r^2^34- Therefore, the total probabihty can be estimated 
as /r^^ which is much smaller than the naive Gaussian estimation L^jr^ 
while the collinear answer L? jr^ is much larger than Gaussian. 

That consideration can be readily generalized for arbitrary number of non- 
coUinear pairs. We expect that F^n oc {Ljr)^'^-^. In accordance with ( |i.5| ) 
the separations Rij should be diminished in the evolution process up to L 
to produce a non-zero contribution into the integral. Suppose that Ryi is 
diminished up to L. Such process (explained in the consideration of the 
pair correlation function) gives the probability {L/ri2f. Next, due to the 
conservation law of the triangular areas all other Rij will increase by the 
factor ri2/L. Then we should diminish, say, i?34 from r-^^ryxji^ down to L. 
Such process gives the probability I? j {ryir-i/^ . Due to the conservation law 
of the triangular areas other Rij will be larger than their initial values by the 
factor ru/ L at the moment. Repeating the process we come to the factor 
(L^/r^)^"^^^ for the n-th order correlation function. Therefore 

A2n = 4n-2 (1.11) 

The above analysis can be generalized for arbitrary geometry. Suppose that 
among the separations Vij are parallel vectors (more precisely, with angles less 
than L^/r^). Let us divide into sets consisting of the pairs of points with 
parallel separations Vij. All points of such set behave as a single separation 
at the Lagrangian evolution. Therefore instead of n we should substitute 
into ( |1.11|) the (minimal) number of the sets. The estimates, obtained above 
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will be supported by rigorous calculations in Sec. |T3. 

Unfortunately, not much can be argued qualitatively about the scaling at 
d > 2. The crucial point for our considerations in d = 2 was the conservation 
of the area. It allowed us to get the correct answers even without calculations. 
In other terms, it is related to the fact that there is a single Lyapunov expo- 
nent at two dimensions. When d > 2 we have only the conservation of the 
d-dimensional volumes and hence more freedom in the dynamics. Consider, 
for instance, the three-point correlation function for non-coUinear geometry. 
Unlike d = 2 we can not assert that it is zero, since now the area of the 
triangle is not fixed and can change during the evolution. Nevertheless, the 
anti-correlation between different Lagrangian trajectories exists, and there- 
fore the answer for the exponent A3 should be larger than 2d which is the 
estimate one would get without the anti-correlation. In the following sections 
we find that = d + {d — l)\/d/ {d — 2). This is determined by the hierar- 
chy of Lyapunov exponents giving the stretching rates at different directions 
— Sec. |iV C| . Note that in the limit of large d the anti-correlation should 
disappear and the answer tends to 2d. The four-point correlation function is 
also determined by a joint evolution of two distances and A4 = A3. 

II. ANALYTIC CALCULATIONS 

We do all the calculations assuming the strain to be delta-correlated in 
time 
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\<7a|3{t^)(7^ily{h)) = D[{d+ l)Sa^Spj, - 5aiy5f3^ - SapS^iy] 5{ti - t2) . (2.1) 



The tensorial structure in is due to isotropy and the incompressibility 



condition diw = tro" = 0. Zero correlation time of the strain allows one to 
derive closed equations for the correlation functions of the scalar!: 

DCF2n{rk) = - E X (k^ - ^jl) i^2n-2(t, Vk') , (2.2) 

DCF2n+i{rk) = - E X (k^ - rj\) F^n-iit, rfcO (2.3) 
- E X3 (In - \r, - r^l, - r^\) F2n-2{t,rk") , 

ijm 

where Vk' is the set with and Vj excluded and r^" is the set with r^, 
Tj and Tm excluded. The dimensionless operator C is written as follows 

■d+ 1 



>^: = E 



ij ij ij 



_,Vrvf, (2.4) 



Eqs. ( PT2| , P73| ) are rather complicated partial differential equations. We 
start our analysis from the pair correlation function. 



A. Pair correlation function 



Due to isotropy and translational invariance, (\Z.'4 ) for the pair correlation 
function can be written as 

^-^-^r'-% {t'^%F,) = -x{r) , (2.5) 

One can easily find a solution of this equation, satisfying the correct boundary 
conditions 
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2 



2 roo dx . , . 



dy xiyjy''-' + f — xiy) 



[2.6) 



d{d- l)D L ^0 ^^v^^^ Jr y 
At r L, the function x(^) is assumed to decay fast enough (say, expo- 
nentiaUy), and it is possible to neglect the terms related to the tail of x(^) 
so that 

Here 

x = ^yx{y)y ~ 

is proportional to the zeroth Fourier harmonics of x(^)- Estimating x — 
x(0)L^ we get F2 ~ x{^){L/rY. It is important that x exists and is non- 
zero, otherwise the answer is different. 

B. CoUinear Geometry 

Here we consider a 2n-th order correlation function of the passive scalar 
regarding that all points ri, . . . , r2n He on the same line. Then as follows 
from during the evolution will remain on a line. The direction of the 



line can be characterized by a random unit vector m{t) with the statistics 
determined by the equation 



following from For the collinear geometry 
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Rijir) = \ri - rj\ exp lj°dt((t) 



[2.9) 



The statistics of the field ( is determined by ( PTT1 , PTB|) which leads toiJi 



{({tiKih)) = Did - 1)S{U - h) . 



[2.10) 
[2.11) 



Using the expressions ( |2T9| , |^.iU| , |^.f i|) we can obtain the closed equation for 
the function F2n{t,rk)' 



(d-1) 



D 



dr 



F, 



2n 



= - T.x{\rz-rj\)F2n-2{t,rk') 



2.12) 



where Vk' is the set with and Vj excluded. 
Let us parameterize the points like 

= ri + e^nli , 



;2.i3) 



where n is a unit vector and k are some coefficients. Then, the equation 
( [k!.i2|) can be rewritten: 

^^^D {dd^ + dl) F2n (e^r,) = - Ex (e^k. - F2„_2 {e^rt') . (2.14) 

This is an ordinary differential equation which has to be solved with the 
following boundary conditions: F2n {e^ri) tends to zero if ^ ^ +(X) and 
remains finite if ^ ^ — oo. The solution is 



F2n{r^) = 



X 



/• + 00 



Ex (e^ln - F2n-2 {e^rk' 



[2.15) 
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If the separations e^\ri — rj \ are much larger than L then the right-hand side 
of (^TT^) can be neglected and we conclude that F2n oc exp(— d^). Thus we 
deal with an extremely strong intermittency when the scaling exponents are 
independent of n. If all separations are of the same order r then we get from 
( [k!.i4|) an estimate 



2n 



(L/rY 



■2.16) 



Note that if the distances strongly differ then it follows from (|27T5]) that it is 
the largest distance that gives the main contribution into 

The analogous procedure can be applied to the odd correlation functions 
of the passive scalar 6. The only difference is that now we should take into 
account also the third-order correlation function of the pumping. Then we 



get 



id-l] 



D 



1-3 , 



F, 



2n+l 



Ex(k* -rj\)F2n^i{t,rk') 



[2.17) 



E X3 {\ri - r-jl, \ri - r^|, \rj - r^l) F2n~2{t,rk") , 

ijm 



where r^/ is the set with and rj excluded and rk" is the set with r^, 
Vj and Vm excluded. Considering all the separations of the order of r we get 
from 



F, 



2n+l 



D [d 



d 

) , 



;2.i8) 



where P3 = X3(0,0,0). The same r-dependence of the odd correlation func- 
tions as in is accounted by the the same structure of the the differential 
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operator in the left-hand sides of i \Z.i'4) and (^TT7|). 



III. DIMENSIONALITY TWO 

As we mentioned above, the 2d case needs a separate consideration because 
of an additional degeneracy of equations for the correlation functions of the 
passive scalar. The degeneracy is associated with the area conservation law 
of any triangle, vertices of which move along Lagrangian trajectories. 

A. Triple correlation function 

As explained in Sec. |T|, the three-point correlation function has a sharp 
peak for the coUinear geometry, whereas for the general position of points 
the answer is determined by the tails of the pumping function and is non- 
universal. Therefore, only the collinear answer is of interest which has been 
already obtained in Sec. |11 B| . Here, we just re-derive the result in a system- 
atic way, starting directly from the equation DCF^ = — Xs- Introducing the 
variables 

a:i = — cosf?, a:2 = — sinf?, s = ri2ri3 sin , (3-1) 
ri2 ri2 



the operator jC ( |2.4| ) can be recast to the following simple formtllll3 

jC = 2x\{d\ + dl) . (3.2) 

Here d is the angle between ryi and ri3, and s is the doubled area of the 
triangle, with vertices in ri, r2, and r-^. Thus, the solution can be easily 
found00 (see also (|4.5|) ) 
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1-00 

-oo 



-oo 



X2 



In 



[Xi - x[y + {X2 + ^2)^ 

^xi - + {x2 - 3:2)2 



X3(ri2,r^23)^3l] 



(3.3) 



One should substitute into ( |3T3| ) the transformations inverse to ( PTT]) 



\ 2^2 



X2 



?^23 



s([a:i — 1]2 + 2:2) 



^2 



(3.4) 



and the analogous relations between r'^2? ^235 ^31 s, x'j^, X2. 

One can easily check, that if s ^ L^, there is no such values of x[ and X2, 
that r'^25 ^135 Q'^d r23 are smaller than L. Therefore, in this case the value 
of F3 will be determined by non-universal behavior of the function ^3 for 
the values of its arguments larger than L. Consequently, properties of the 
correlation function F3 are non-universal, in agreement with the qualitative 
discussion of Sec. |. 

Let us consider F3 at s ^ L^. Since we assumed that both ru and ri3 
are much larger than L, the condition s ^ gives the inequality 7? < 
F^/ (^12^13) 1, that is we consider geometry close to collinear. In this 
case, the main contribution into the integral ( |3.3| ) is made by the region of 
integration, where all r' are smaller than L. In particular, r^2 ~ ^5 which 
allows to estimate X2 ^ s/L^, which is the same as X2/X2 ^ ^12/^^ ^ 1- 
Therefore, we can expand the resolvent in ( |3.3|) and write 

X2 fdx2 fX3{r'i2^r[^,r'2-i)dx[ 



(3.5) 



27tD J x'2 ^ (xi - x'l)^ + x'i 
In the limit s <C L^, it can be further simplified since the main contribution 
to the integral ( |3.5| ) is associated with the region x'2 <C min(l, x\): 
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1 foo B^L'^ 

^3 = 7^ /n XsiruL ri3^, r230 « 777^ — 2 — 2^ • (3-6) 

DJO 2L»max (rf2,rf3,r^3) 



The expression ( p.t)|) is in accordance with the estimate (|27T§). Note that 



( p.6|) has no singularity when any two points coincide as long as at least one 



distance is finite. 

B. Four-point correlation function 

In this subsection we derive the result for the four-point correlation func- 
tion starting directly from (\2.'4 )- Again, there are two regimes for which one 



can find the answer. For the collinear geometry, the consideration is very 
similar to one done in the previous section and reproduces the result 



Here we will find the answer for the non-coUinear geometry. Note that its 
estimate is already known from Sec. |. The equation i\Z.'Z^ for the four-point 
correlation function F4 is 

—DjCF^ = x(^i2)^2(^34) + permutations . 

The property of the operator (characteristic of the large-scale advect- 
ing velocity) is that the solution of this equation is reducible into pieces, 
corresponding to each term in it's right-hand side: 

F4 = F4(ri2, r34) + ^4(^34, ^^12) + ^4(^^13, ^^24) 

+F4(r24, ris) + F4(ri4, r23) + ^4(^^23, ru) , (3.7) 

To find F4 we should solve the equation 
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(3.8) 



In terms of the variables (f? is the angle between ru and vm) 

xi = — cost;', X2 = — smv, s = ri2r34 sm , 
ri2 ri2 



the operator C for d = 2 i 



(3.9) 



£ = 2x^(9^ + 92^). 
The solution of ( fTBI) can be written as a double integral 



(3.10) 



1 /•+00 

SttD J-oo 



> . 



dx2 

X2 



x{A2)F2{r'M) In 



{Xi - x\Y + {X2 + x'2) 



I \2 



xi — x'lY + (3:2 — x'2 



i\2 



ri2 = 



\X'2 



^34 — 



\ 



s{xi + X2) 



x'2 



(3.11) 



We shall calculate the integral in the limit s ^ when there are several 
simplifications. First, since r'12 ~ L, we can write X2 ^ s/L^. Hence, r'^^ > 



sx'2 ^ s/L ^ and we can use the asymptotic form ( P77|) of F2. Second 



like for the three-point correlation function one can show that X2 ^ 2^2, and 
one can expand the logarithm. Finally, we can write 

00 dx'i 



X2X2 



-00 



\ ^2 / °° Wl + ^2 ] [{x'l ~ 3:1)2 + X2] 



47r2L)2s ^0 

Integral over x'l can be easily calculated and we get 

-00 dx'o 



(3.12) 



X2X2 



1 



2tiD^sJ^ x'2 xl + Ax'i^\\x'2, 
If xi is not anomalously large, we can disregard it in the integrand, and find 



(3.13) 



Fa = 



X2CX2 



c = j^ x(Ord^ 



(3.14) 



IV. DIMENSIONALITIES LARGER THAN TWO 



Here we treat correlation functions of the passive scalar foi d > 2. In this 
case, the degeneracy inherent to = 2 is absent, and the consideration is the 
same for all which is thus considered as a parameter. Of course, direct 
physical meaning can be attributed only to d = 3. 

We will calculate the three and four-point correlation functions. Exactly 
as it was for d = 2, the operator jC has the same form for both correlation 
functions. Namely, we have to solve the following equations 

-m = X3, -^F, = x{n2)F2{rM) , (4.1) 



where F4 is defined by (|X7|) . Then, we can introduce the variables (|371| ) for 



F3 and (|3.y| ) for F4. In these variables the operator C has the following rather 



simple form 

C = dxl{dl + dl) + (d - 2){dl + ddr) , r = In (s/L 



Therefore, in order to solve Eqs. ( |4.i| ) we have to find the resolvent IZ of the 



operator C which satisfies the equation 

—CTZ = 5{xi — Xi)5{x2 — X2)5{t — r) . 

and the following boundary conditions: First, TZ should go to zero when 
xi ±(X), X2 +CXD, X2 0, and r +00. Then TZ should tend to 
a constant at r ^ — (X). It is more convenient to work with Hermitian 
operators, therefore it is useful to make a substitution 
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X2 

IZ = ——r exp 

ax2 



d 

--^r-r 



R{xi, x[, X2, X2, r, t') 



Then we obtain an equation 

^d^R d'^{x2R) d-2fd^R 



X 



^ dxf 



+ X2- 



dxo 



+ 



d \dT^ 4 



R 



= —6{xi — x'i)6{x2 — X2)6{r — t'] 



It is natural to seek the solution in the following form 



R = 



1 



dk 



da 



— I — a, a:2, ^2)6 



/ \ ik{xi—x'i)+ia{T—T') 



(4.2) 



The function w satisfies the following equation 

nd'^u du 



X 



dx2 



{k^x\ + v'^)u = —\[j^2^{^2 — x'2) , 



where 



V = 



(d-l)2 d-2 



+ 



(4.3) 



The solution of this equation can be readily expressed in terms of the Bessel 
functions of imaginary argument (see e.g. El) 

w = {e{x2 - x'2)K,{\k\x2)h{\k\x'2) + ^(4 - X2)K,{\k\x'2)h{\k\x2)} 



Here 9{x) is the step function, equal to one if a; > and zero otherwise. 

Now we should substitute u back into Integral over k can be cal- 

culated analytically with the help of Eqs. (6.672, 8.820) fromlil. Then we 

get 
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R = 



(iar(z/ + 1/2) 



2^x2X2 27r V{v + 1 



xF 



V 3 1/ 1 
- + -,- + -; + 1; 
2 4' 2 4' ' 



2 + 
2x2^^2 



X2 I ^2 I V"^! — "''I 
2 



r\2 



X2 + 2:2^ + (3:1 — x'l)'^ 
Here F{a, P;j;x) is the hypergeometric function and T{x) is the Euler 
Gamma function. 

In d = 2, the integral over a is trivial, and the resolvent can be easily 
reproduced^'li^ 



exp [ia[r — r 



(4.4) 



7^ = 



1 



Sirx'o^ 



In 



Xi — x'l) + {X2 — x'2 



bir - t') 



(4.5) 



[Xi - Xi)2 + {X2 + X'2)^ _ 

Convolution of the resolvent with the right-hand sides of Eqs. ( ^TTp depends 
on the properties of the sources. Therefore, below we consider separately 
both correlation functions. 



A. Three-point correlation function 

The solution of ( |4.i|) for the three-point correlation function can be written 
in the following form 

^3 - r d4 r d^'i r (^1. <^ ^. ^'^"'''''y^-'^ . (4.6) 

J— 00 J— 00 ) 



The relations between the variables ^125^135^23 t' ^rj' are as in ( |a.4j) . 
Remind, that r = In s/L^. 

Like it was in d = 2, the behavior of F3 is very different for the cases of 
7? <C 1 and 7? ~ 1. Let us first consider the case of not very small angles, 
namely 'd ^ i^^/(ri2ri3). Since both ri2 and ri3 are much larger than L, the 
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area s of the triangle is much larger than L^, which means that r ^ 1. On 
the other hand, since xa decreases very rapidly when any of rij is larger than 
L, the area s' can not be much larger than L^. Therefore, r' is of the order 
unity in the integral ( [4.6|) . Thus, we see that r — r' is always positive and 



much larger than unity. On the other hand, from the condition i!) ^ L'^/r'^, 
it is easy to check that for a typical configuration, contributing into (|^ ) the 



condition InA <ti r holds. Therefore, shifting the contour of integration in 
( [4.4|) into the upper half-plane, we will meet a branch point of the integrand. 



that originates from u and is situated at a = i{d — 1) Y<i/ {d — 2). Because of 
the large value of r, the integral will be determined by a near vicinity of the 
branch point. Therefore, we can write 

A3 



71 oc exp 



-T 



(4.7) 



A3 = d+{d-l)^d/{d-2). (4.8) 
Substituting the expression into (|4T^ ) we get 



F3 ^ iPs/D) {L/rf' , (4.9) 



Note that A3 > 2d. 

Let us now consider the limit i!) <C L'^/r'^ and reproduce the coUinear 
result. In this case, as we shall see, the main contribution to (|4?^ ) is made 



by <C 1. The resolvent IZ in the limit t?, t?' <C 1 can be found directly from 
the representation (|4T2|). The smallness of the angles implies that X2,X2 <C 1. 



Using the asymptotic expansion of the Bessel functions one gets 

22 



u = 



2v 



1 (4,0) 



We see, that in the main approximation the /c-dependence disappears of u 
and we can integrate over k in ( |4.2| ). Then we get 



R 



5{xi — x'l) roo da 



exp 



laiT — T 



2^x2X2 2'Kv{a] 
The integral ( ^.iil) can be calculated analytically 

d 



v{a) 



In — 

X2 



(4.11) 



R 
X 



1 



X2 



7r^x|d(d-2 
d- 1 



-Kq {X) exp 



(r - r') 



^ ( /^2 I 1 2 

2 V"^^"""^^ 4 



(4.12) 
(4.13) 



The (5-function forces the ratio ^22/^13 to be equal to ryil^xz- Integration 
over one of the distances, say makes both of r' to be of the order L (we 
believe that ryi ~ ^13)- Let us consider the integral over the angle -d' . It is 
easy to see, that the argument of i^o in (|4.i2|) is always large. Therefore we 
can use the asymptotic form of this function and write 

d 



^1/2 



In 



j2 ^ 



d- 1 



In^ 



+ 



d 



2 \ d-2 



ly 2 



The main contribution to the integral is made by the vicinity of = dr^ jU 



Hence F3 



L 



From the assumption f?' <C 1, we see that there should 



be 7? <C L? jr^^ otherwise the main contribution comes from 7?' ~ 1, and the 
expression for the resolvent ( [4.12| ) is inapplicable. 
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B. Four-point Correlation Function 



From ( |4. i|) it follows that the answer for the four-point correlation function 



can be written in the form (|5.V| ) where 



~ 1 fOO f + OO f + OO 

^^4(^-12,^34) = L / / dr 

) JO J— 00 J— 00 



x7^(xl, xi, X2, X2, r, r ) x(^i2)^2(r34) , (4.14) 

The variables r'12 and are expressed via , X2 5 and r' by ( [TTTl) . The anal- 
ogous relations hold for the variables ri2 and r^^. It can be more convenient 
for the present purposes to pass from the integration over x[, x'2^ and t' to 
the integration over 7?', r']^2) ^34- Then ( |4.14|) can be replaced by 



The expression for A can be written as follows 

2 sin sin / , . 
— M X5) 

^^12^^34/ (^12^34) + ri2r34/ (ri2r^4) - 2 cos 7? cos 'd' ' 
The case i!) <C L'^/r'^ can be analyzed in a way, presented in the previous 
subsection for the three-point correlation function, leading to the expression 
( [k!.lb|) and to the law (|27T^). Below we analyze the case i!) ^ L'^/r'^. In this 



case, there are two different regions in the integral over ^2 in the integral 
making contribution to F4. The first region is determined by the inequality 
^2 ~ 1, which corresponds to r'^^ ~ L. By the order of magnitude, this 
contribution is equal to the three-point correlation function, because in this 
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case, the value of t' in ( [4. i4|) is of the order of unity, while r 1 and 



therefore all the arguments presented in the previous subsection are valid. 
The second region is defined by the condition ~ 1- It corresponds to large 
values of r34, for which we can use the asymptotic behavior of F2, given by 
( ^.7|) . Another simplification is that in this region r'^^ ^ therefore 



we can believe 

A ^2smdsmd' ^-^^ < 1 (4.16) 
^12^34 

Then, we can put the hypergeometric function to be equal to 1, and write 



X 



U V 

The integral over angle d' can be easily expressed via the Euler F-functions. 

The integral over v should be calculated with the cut-off on L. Therefore, 

we get the following integral over a 

^ r(z/+l/2)r(z//2-l/4-za/2 + d/4) /ri2r34 . .V""'^^' 

F4 oc / da^- , ' - , rrr^ — sm if 

J r(z/+l) r(z//2 + l/4-za/2 + d/4) V I 



d/2 + v +l/2 + ia-l u^^ ' \L^, 
Again, we can shift the contour of integration up to the branch point de- 
termined by the same condition z/ = 0, and we get the same answer as for 
the three-point correlation function. Thus, both contributions possess an 
identical r-dependence giving 



P*--F^A-) . (4-18) 
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with the same exponent (|4.^|) . 

The contribution ( [4 .181 ) to F4 is the leading one only if d > \f2 + 1. If 



d < + 1, then along with (|4.i^ ) there appears an additional contribution 
into F4 due to a pole of the integrand in ( [4.iY|) . It originates from the zero 
of the denominator d/2 + z/+l/2 + iQ;, existing only at d < \J2 + 1. The 
contribution behaves like 

, j^.d{d+i)/{d-i) 
oc - . (4.19) 



r 



Comparing the law ( |4. with ( [4.i9| ), we conclude that in the region of its 



existence, that is at d < 1 + a/2, the term ( [4. is the leading contribution 



to F4(ri2, T34). Particularly, this is the case in d = 2, where the contribution 
( |4. i9|) behaves as (L/r)^, in accordance with (|3.i4| ). 



C. Instanton 

To understand better the underlying Lagrangian dynamics let us outline 
briefly another way of calculation, based on the fact that they are rare events 
that contributes the correlation functions at large scales. Therefore, some 
kind of instanton formalism can be applied, the main task here is to recognize 
the relevant degrees of freedom. In this way, we shall establish a relation 
between the scaling exponent A3 and the Lyapunov exponents of the smooth 
flow. The Lagrangian distances Rij are all determined by a single matrix 
W = T exp{J adt) via Rij = Wvij. To find the correlation functions, we 
should be able to average over the statistics of the matrix W. The way to 
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do it was proposed in^, we shall followlll and write 



W = RT, 



(4.20) 



where R is an orthogonal matrix and T is an upper-triangular matrix: Tij = 
if z > j. The matrix R can be excluded from the consideration due to isotropy, 
then, representing T as 



T,i = exp(pi), T,j = exp{pi)r]^j if i < j 



(4.21) 



we can write the action describing the stochastic dynamics of p and ij 



i=l 



^tp^ + D 



d{d -2i + r 



iD 

+ T 



+iDd exp(2p, - 2pj)/i^^- + 2iDd I^ijf^ik exp(2pfc - 2pi)r]kj 

i<j i<k<j 

+ Y f^tjdtmj + 'iDd Y IJ'imlJ^inVkmVkn exp{2pk - 2pi) 

i<j i<k<m,n 



(4.22) 



Here rrii and Pij are auxiliary fields, conjugated to pi and r]ij respectively. 

The variables pi describe stretching of volume elements in the flow, while 
ijij describe the direction fluctuations of a given vector with respect to the 
main axis of the strain matrix a. Note that the constants = d{d — 2i + l)/2 
entering the action ( [4.22|) are the Lyapunov exponents. Using (|4.22|) , we can 
rewrite F3 as follows 



F3 = / V pVrjV pVm exp i ( dt C + In { [ dtxz 

J J—oo \j—oo I 



(4.23) 



The variables r]ij are irrelevant for the evaluation of the scaling and are 
only responsible for the angular behavior. Since all three Rij always lie on a 
plane, we can leave only two of pi and write an effective Lagrangian 
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iD 



C = ruaidtpa + Afl) + mb{dtpb + A^) + — [d{ml + ml) - {ma + mt 



Then, the dependence on L/r can be extracted from the following expression 



J Vpa^bVma^b exp I J_ dtjC + In { f dt xs 



-oo 



■oo 



(4.24) 



Now we can write the instanton equations for the extremum of the exponent 
in i WTI^ 



dtPa + Aa = -iD[{d - l)ma - mb] , dtpb + A& = -iD[{d - l)mb - ma] , 



F a Pa F dpb ^-oo 

The boundary conditions are ma^b ^ as t 
Note that the "energy" 

D 



dtx3 (4.25) 
— oo and pa^b = at t = 0. 



E = limaK + rubh 



[d{ml + 



2^ 



JUa + mb) + —X3 

^3 



1 



(4.26) 



is a constant. From the boundary conditions we deduce that it should be 
zero. 

Let us explain the qualitative behavior of the solution. We consider the 
evolution backwards in time. At small times all Rij are large, and therefore 
the derivatives of xs the right-hand sides of Eqs. ( [4.25|) are zero. Hence 
ma^b are constants such that pa^b diminish and Rij also become smaller. Then, 
at some moment all Rij become of the order of L. Then, derivatives of xs 
can not be disregarded and during some short interval of time when Rij ~ L 
the values of ma^b will change to vacuum zero values. The derivatives of pa^b 
change sign, and Rij start to grow. Note, that if only one of Rij is of the order 
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L and the other are still much larger than L, then ^3 is small, the derivatives 
in Eqs. ( [4.2b|) are ineffective, and the solution will never reach its vacuum 



stage. Thus, we should tune the conditions so that all Rij will become of 
the order L simultaneously. Finally, we come to the set of conditions for the 
initial stage 

dtPa = A = -Xa - iD[{d - l)ma - mfe] , 



dtPb = A = -\b - iD[{d - l)mb - 

L 
'2 



E = lipiaXa + rnbXb) - ^ [(i(mj + ml) - {ma + m&)^] = 



From here we find the value of A 

1 



A = 



(Aa + A,)2 + ^(A«-A6)2 (4.27) 



2V "" '^ 

Calculating the action, we find that F3 ~ (L/r)^3 with 



Af = ^- — — (4.28) 

The value of A3'^ is minimal (that is F3 is maximal) if we take the two 
largest Lyapunov exponents: A^ = Ai = d{d—l)/2 and A5 = A2 = d{d — 3)/2. 
Substituting it into ( |4.:^^ ) we reproduce ( ^TB|) . The instanton found has a 



long life-time, proportional to ln(r/L), therefore the above consideration has 
to be valid also for a velocity finite-correlated in time. 



D. Discussion 



We thus see that many-point correlation functions are not scale invariant 
because of strong angular dependence. One may consider averaging over 
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different geometries, for instance, integrating over the angle between any 
two vectors keeping R"^ = J2 rfj in 2d. As a result of such averaging, the 
object appears which depends on R only and is thus scale invariant. Does 
such averaging restores also the normal scaling? One may notice that the 
main contribution into the angular integral gives the region of small angles 
near collinear peak i!) ^ (L/r)'^; since there are n — 1 angles in the n-point 
correlation function then one gets F„ oc r~^" that is normal scaling is restored 
indeed. That means that the increase at small angles and decrease at large 
ones (relative to a normal scaling) are of the same order and both caused 
by the same mechanism. It is unclear what is the way - if any - of natural 
average over geometries that restores normal scaling in d > 2. 

We thus discovered an intermittency build-up in the direction opposite to 
the cascade. It is instructive to compare this with an intermittency discovered 
at small scales when the cascade flows upscales in a compressible flowEl. 

What will be for a finite correlation of a in time? It is clear from Sects. 
0,0 that the dependence r^^ both for the pair correlation function and for 
the correlation function of any order for collinear geometry is valid as long 
as the correlation time of a is much less than In {r / L). Under the same 
assumption, all the results obtained in two dimensions will be valid (up to 
a single numerical factor in front of any correlation function) for a finite- 
correlated strain as well. As far as higher dimensions are concerned, it is 
clear that some anti-correlation between contraction of different distances 
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will be present, it is likely that it will be governed by the same exponent 
A3. Indeed, ( |4.28| ) has to hold in a finite-correlated case as well, A3 is 



thus determined by two largest Lyapunov exponents which are likely to be 
proportional to {d — 1) and {d — 3) respectively. 



V. NON-SMOOTH VELOCITY FIELD 



The advection of the passive scalar by the non-smooth velocity in the 
framework of the Kraichnan model is described by Gaussian velocity field 
with the pair correlation function 



{t, r)w'^(0, 0)) = 6{t) {"l/o(5"^ - r-^[{d + 1 - j)r^6''^ - (2 - 7)r V] } 



Here 7 is a measure of velocity non-smoothness, < 7 < 2. The generaliza- 
tion of the equations ( [2.2| , [2.3|) for nonzero 7 needs replacing jC by 



E r,7[(d + 1 - 7)r?/"'^ - (2 - 7)^^^] V^Vf • 

The pair correlation function is now F2{r) oc r'^~'^i. 

An interpretation of the extra factor comparing to ( |2.7|) is related to 
the fact that every Lagrangian distance R generally grow by a power law t'^^^ 
as distinct from an exponential law at 7 = 0. In other words, stretching is 
uniform in the logarithm of scale at 7 = and decelerating at 7 7^ 0. 

We cannot yet find the high-order correlation functions at arbitrary 7. 
Fortunately, at the limit of extremely irregular velocity 7 = 2 the operator 
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turns into 

c, = -{d-i) [(Ev.)'-Ev?^ 

For translationally invariant functions the term (EV^)^ can be discarded. 
We are thus left with diffusion equation, it is straightforward to show that 
if the pumping is Gaussian then the scalar statistics is Gaussian too (both 
at the scales larger and smaller than the pumping scale). If the pumping 
is non-Gaussian, then the scalar statistics is getting Gaussian at the scales 
distant form the pumping scale (odd moments and cumulants of even mo- 
ments decrease with the growth of r/L faster than the respective Gaussian 
moments). The third moment, for instance, decreases with r/L faster than 
F2^\ this can be shown by substituting the resolvent {d > 2) 

into the expression 

fP = J dxidx2dx3 n{xi, X2, X3, 2/1, 2/2, 2/3) Xsixu, Xl3, X23) 

C3 



[yh + yh + yhV-'' 

C3 = ^'Jd^^t'i)^ I dCdvXs (2I77I, I V3C + Vl I V3C + 77l) (5-1) 
One sees that oc r'^~^^ which is decaying with r much faster than F2^^ oc 

^3-3d/2 

At 7 < 2 one can directly check that the scalar statistics is non- Gaussian 
even for a Gaussian pumping. Employing perturbation theory with respect 
to ^ = 2 — 7, one may try to prove at least that the scaling is normal and 
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the angular anomaly is absent at < 7 < 2. Let us do this for the triple 
correlation function. The operator L to the first order in ^ is £ = £0 + 
with 



5°^[(d-l)lnr., + l]-^ 



Then, we should find t\Ff^ and integrate it with the resolvent. We have 

^ (0) _ _ [ 3(d-l)(3d-2) 
" ~ n[r?2 + r?3 + ri3]^ 



d{d-\) 



[r?2 + r?3 + risl^+i 
2d[d - 1)2 



(^13 - ^23)^ (^12 - ^23)^ (^12 - ^13)^ 



^12 



^13 



/v>2 
^23 



[^12 + ^13 + ^23] 



rTo In 



2 Irf+l '12 



2 2 
^13' 23 

^12 



+ ^13 In 



^12^23 I 2 1 ^12^13 1 



There are three kinds of terms here. From the symmetry reasons, it is enough 
to know the values of the following integrals 



-I 



dxidx2dxs 



2i3d/2-l 



2 — 



3 — 



/ 



(a;?2 + ^13 + ^23)^ [{xi - ViY + {x2 - 2/2)2 + (xs - 2/3)2] 

(xfg — X2'^'^dxidx2dx^ 

AiiAi + ^?3 + A^y-^' [(^1 - 2/1)' + (^2 - 2/2)^ + (^3 - v^?f"^-' 

a:22 In ^dxidx2dxs 



(X?2 + X?3 + xhY^^ [(Xi - 2/i)2 + ix2 - 2/2)^ + (^3 - 2/3)']''/'"' 

With the logarithmic accuracy one finds 

^1,2,3 



^1,2,3 - 



ivh + yh + ^23)'^ 

27,3d/2 



31 In 



Vu + 2/13 + 2/23 



L 



(5.2) 



3((i- l)r(3d/2- 1) 



2^3d/2 



3d{d- l)r(3(i/2- 1) ' 



A. = 



TT 



3d/2 



3d{d-l)T{3d/2-l) 
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The result can be expressed in terms of /i,2,3 



In 



C3r{3d/2 - 1) 
L 



[{9d - 6)/i - 3dl2 + 2d{d - l)h] 




L 



This is the first term of the expansion with respect to ^ of the function 
F3 oc {yf2 + ?/i3 + 2/23)^^^ oc r^^^^^ which is the normal scaling for the triple 
correlation function. 

Analysis of the integrals /i,2,3 shows that there is no angular singularity 
at coUinear geometry. By a direct calculation one can check that the scal- 
ing of the triple correlation function is the same for three points on a line. 
This is natural since nonzero 7 destroys degeneracy, collinearity is no longer 
preserved during the Lagrangian evolution so the correlation functions at 
coUinear geometry have no anomaly, similarly to what has been established 
by 7-expansion at small scalesS'S. 



We have studied the correlation functions of a passive scalar in the frame- 
work of the Kraichnan model on distances larger than the scalar's pumping 
length. In the Batchelor limit, the collinear anomaly has been found: scal- 
ing behavior of many-point correlation functions for the collinear geometry 
(where some points lie on a line) strongly differs from one for general ge- 
ometry. The anomalous scaling is observed in the interval of angles which 



VI. CONCLUSION 
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decreases with increasing scale. This violation of a conventional scaling be- 
havior is related to a strong correlation between different Lagrangian trajec- 
tories occurring in the Batchelor case that is for distances smaller than the 
viscous scale of the velocity field. For larger distances (at the inertial interval 
of turbulence) the scale invariance of scalar statistics (yet not Gaussianity) 
is likely to be restored (remember that we consider the scales larger than the 
scale of scalar 's pumping) as it is confirmed by our calculations in Sect. [V]. 
At even larger scales (beyond velocity correlation scale that is an external 
scale of turbulence) the scalar statistics has to be Gaussian. 
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